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Abstract 

Measurements  have  been  made  of  the  ambient  magnetic  field  and  magnetic 
field  gradients  near  the  Naval  Undersea  Center  (NUC)  Oceanographic  Research 
Tower.  The  experiment  and  Its  results  are  summarized  In  this  report.  The 
NUC  Tower  Is  located  approximately  7/10  of  a mile  off  the  California  coast 
near  San  Diego,  and  Is  the  proposed  site  for  the  shallow  water  trials  of  the 
ARPA  Internal  Wave  Magnetic  Sensing  (IWMS)  experiment.  The  measurements 
described  here  were  made  in  order  to  accurately  determine  the  magnetic  field 
and  associated  gradients  of  the  NUC  tower,  so  that  their  significance  as  a 
possible  source  of  noise  and  Interference  during  the  IWMS  experiment  may  be 
assessed.  An  analytic  model  of  the  magnetic  field  of  the  NUC  tower  which 
accurately  describes  the  field  and  associated  gradients  Is  also  described. 
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I.  Introduction 


The  Internal  Wave  Magnec-C  Sensing  (IWMS)  experiment  will  utilize  an 

1 2 

extremely  sensitive  suj^ erconducting  magnetic  gradiometer.  ’ The  successful 

implementation  of  that  experiment  to  measure  the  magnetic  field  gradients 

3 

generated  by  ocean  waves  requires  the  fielding  of  the  superconducting 
gradiometer  on  a suitable  platform  in  an  ocean  environment.  Above  surface 
observations  of  shallow  water  Internal  waves  have  been  proposed  to  be  con- 
ducted on  the  Oceanographic  Research  Tower  of  the  Naval  Undersea  Center  (NUf ) 
in  San  Diego,  California.  The  NUC  tower  is  located  approximately  7/10  of 
a mile  offshore  and  has  been  the  site  of  several  internal  wave  experiments 
in  the  past.  However,  the  tower  is  constructed  primarily  of  ferrous  materials 
and  can  be  expected  to  significantly  alter  the  earth's  magnetic  field  in  its 
vicinity.  The  superconducting  gradiometer  developed  for  the  IWMS  experiment 

2 

has  very  low  noise  characteristics  in  the  ambient  earth's  magnetic  field, 
but  operation  in  the  presence  of  a disturbed  magnetic  environment,  such  as 
that  near  the  NUC  tower,  poses  a number  of  problems,^  which  must  be  addressed 
for  the  successful  completion  of  the  IWMS  experiment.  In  order  to  accurately 
assess  the  effects  of  the  tower's  magnetic  environment  on  the  IWMS  experiment, 
measurements  were  made  to  determine  the  magnetic  field  and  magnetic  field 
gradients  in  the  vicinity  of  the  NUC  tower. 

This  report  sunmiarlzes  the  concept  and  execution  of  the  experiment  and 
presents  the  results  for  the  measured  magnetic  field  and  associated  gradients. 
Those  results  are  then  used  to  develop  a simple  analytic  model  of  the  magnetic 
field  near  the  NUC  tower.  This  model  is  of  sufficient  accuracy  to  use  with 
confidence  in  the  calculation  of  magnetic  field  properties  in,  and  near,  the 
region  where  the  experiment  was  conducted. 
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II.  Experiment  Summary 

The  objective  of  the  experiment  was  to  measure  the  ccmponents  of  the 
steady-state  magnetic  field  and  Its  gradient  matrix  at  selected  points  near 
the  NUC  ocear ographlc  tower.  Nine  points  on  the  west  side  of  the  tower 
were  chosen  which  span  the  region  where  the  boom-mounted  superconducting 
gradlometer  would  be  located  during  the  shallow  water  trials  of  the  IWMS 
experiment.  The  locations  of  these  points  are  Indicated  schematically  In 
figure  1.  To  obtain  the  field  Information  required,  a three  axis  fluxgate 
magnetometer  was  used.  This  directly  measures  the  three  orthogonal  com- 
ponents of  the  magnetic  field.  A single  axis  fluxgate  gradlometer  having 
a longitudinal  configuration  (coaxial  pick-up  loops)  was  used  to  measure 
the  magnetic  field  gradients.  As  will  be  discussed  below.  In  order  to 
j measure  all  the  Independent  components  of  the  gradient  matrix,  data  must 

be  taken  for  a number  of  angular  positions  of  this  single  axis  gradlometer. 

A special  turntable  mount  was  designed  and  built  which  allows  the  gradlometer 
to  be  remotely  rotated  about  two  axes  and  locked  In  position  at  any  of 
several  angles.  Figures  2,  3 and  4 show  the  complete  magnetometer,  gradlo- 
meter and  mount  deployed  at  one  of  the  positions  during  the  course  of 
the  experiment. 

Prior  to  the  experiment  on  the  NUC  tower,  the  magnetometer  and  gradl- 
ometer were  calibrated  and  balanced  at  a field  site  at  the  La  Posta 
Astrogeophyslcal  Observatory,  located  approximately  70  miles  east  of 
San  Diego.  The  calibration  coefficients  determined  there  were  used  through- 
out the  tower  experiment.  The  gradlometer  balance  was  not  stable  from  day 
to  day,  so  It  was  redetermined  for  each  measurement  point  In  the  data 
analysis  (this  Is  described  In  Section  III).  The  experiment  procedure  for 


2 


obtaining  the  required  data  (described  in  Section  II-2)  was  also  carried 
out  at  the  La  Posta  field  site  using  a knovm  magnetic  field  and  field 
gradient,  in  order  to  check  recording  and  analysis  techniques. 
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Figure  1.  Geometry  of  NUC  Tower  Magnetic  Field  and  Field  Gradient  Measurements 


Fig.  2 Fluxgate  magnetometer,  gradiometer , and  turntable  mount  on  a fiberglass  boom 

deployed  to  the  West  of  the  NUC  Oceanographic  Research  Tower.  The  instruments 
are  shown  at  the  C-3  position  indicated  in  figure  1. 


Close-up  view  from  the  water  of  the  instrument  package  deployed  at  the  C-3 
position.  The  turntable  is  positioned  at  0°  tilt  angle  and  0°  rotation  angle 


Fig.  4 Close-up  view  of  instrument  package  from  the  upper  dock  of  the 
NUC  Tower.  The  turntable  is  positioned  at  45”  tilt  angle  and 
90”  rotation  angle  in  this  view. 
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II-l  Magnetic  Field  Gradient-Measurement  Rationale 

The  gradient  of  the  magnetic  field  may  be  written  as  a second-rank 
tensor,  G.  In  dyadic  notation 


G 


7B 


(1) 


where  B Is  the  magnetic  field. 

In  a particular  coordinate  system,  which  Is  defined  by  the  set  of 
orthogonal  units  vectors  where  1>1,2  or  3,  the  magnetic  field  may  be 
represented  as  a 3 component  vector  and  Its  associated  gradient  as  a 3x3 
matrix.  The  Individual  elements  of  the  gradient  matrix  are: 
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x^*G*Xj 


(Xj^*7)(B‘Xj) 


(2) 


Bj  Is  the  component  of  the  magnetic  field,  B,  along  the  x^  direction. 

From  the  time- Independent  Maxwell's  equations  In  free  space,  both 
the  divergence  and  curl  of  B are  zero. 


7*B  - 0 , 

(3) 

7xB  - 0 . 

(4) 

These  require  that  G be  traceless  and  symmetric,  respectively.  Consequently, 
for  these  conditions  there  are  only  5 Independent  components  of  the  (3  x 3) 
gradient  matrix,  G . In  order  to  determine  all  five  Independent  components 
using  a single-axis  gradlometer,  measurements  must  be  made  at  a number  of 
different  angular  orientations.  The  transformation  law  for  rotations  of  the 
gradient  tensor  may  then  be  used  to  evaluate  all  the  Independent  components. 
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Specifically,  the  gradient  tensor,  G,  transforms  under  a rotation. 


described  by  an  operator  R,  according  to 

G'  - R G r"^  . (5) 

SS  1 

Any  proper  (physical)  rotation  is  equivalent  to  three  successive,  independent 
rotations  such  ae  those  defined  by  the  Eulerlan  angles.  Consequently  the 
rotation  operator  R may  be  written  as  the  product  of  three  separate  operators, 
each  only  dependent  on  one  of  the  Eulerlan  angles.  However,  for  a gradlometer 
in  a longitudinal  configuration,  only  two  Independent  angles  are  required  to 
specify  all  the  physically  distinct  orientations  of  the  gradlometer.  (The 
third  independent  angle  may  be  taken  as  corresponding  to  a rotation  about 
the  symmetry  axis,  and  the  output  from  an  ideal  axial  gradlometer  is  invariant 
under  such  a rotation.)  Therefore,  in  the  following  discussion,  only  two 
independent  rotations  will  be  considered. 

In  order  to  provide  a simple  physical  representation  basis  for  the 
magnetic  field,  associated  gradients,  and  rotation  operators,  a fixed  coordinate 
system  defined  with  respect  to  the  NUC  tower  is  used.  This  coordinate  system 
is  indicated  in  figure  5.  The  axis  is  pointing  downward,  along  the  center- 
line  of  the  tower,  the  axis  points  East,  and  the  x^  axis  points  North. 

The  origin  is  taken  to  be  on  the  tower  center-line,  at  the  ocean  floor. 

The  turntable  on  which  the  gradlometer  and  magnetometer  were  mounted 
is  capable  of  rotation  about  2 angles,  which  are  also  defined  in  figure  5. 

The  tilt  angle,  0,  corresponds  to  a rotation  about  the  x^  axis  in  a positive 
sense,  and  defines  an  Intermediate  coordinate  system,  with  axes  labelled 
as  Xj"  , Xj"  , and  x^"  (X2"  “Xj).  The  rotation  angle,  i(<,  corresponds  to  a 
rotation  about  the  x^"  axis,  but  In  a negative  sense.  The  combination  of 
these  two  rotations  then  defines  a coordinate  system  labelled  as  x^'  , x^' 
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and  x^'  (x^'  “ x.^' ) which  is  fixed  to  the  gradiometer  and  magnetometers, 

The  corresponding  rotation  operators,  which  may  be  expressed  as  3 x 3 
matrices  in  this  representation,  are 

cos  0 0 

0 1 


and 


0 


sin 

0 

0 

cos 

-sin 

sin 

cos 

0 

0 

-sin  0 
0 

cos  0 


(6) 


(7) 


where 


R * R j R^ 

4i  0 


(8) 


The  gradiometer  was  initially  aligned  (0*0  and  » 0)  such  that  it 
measured  the  0^2  component  of  the  gradient  matrix.  Consequently,  it 
measures  the  component  when  rotated  through  the  angles  0 and  4». 
Evaluating  the  component  of  the  gradient  tensor  from  equation  (5), 
one  finds  that  Che  result  may  be  written  as: 


• X*  • G*  X* 

22  *2  X 2 


“ 0^2  II '•"I  cos  + Y cos  20  (1  - cos  2i|<)] 


+ Gj^j^  cof.  20  (1  - cos  2\l>)  - Gj^2  sin  20(1  - cos  2\i>) 


+ (Gj^2  COS0-G2J  sin0]  sin  2i|; 


(9) 
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Equations  (3)  and  (4)  have  been  used  In  order  to  reduce  the  number  of 


gradient  matrix  elements  appearing  in  (9)  to  only  five  independent  ones. 
It  is  also  clear  from  this  result  that  all  five  independent  elements  are 
required  to  completely  specify  the  observable  arbitrary  angles  0 

and  \i>.  Consequently,  measurements  made  at  a suitable  number  of  angular 
orientations  allows  the  determination  of  the  complete  gradient  matrix  in 
a fixed  coordinate  system. 


I 

i 
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Figure  5.  Coordinate  systems  and  angles  describing  the  turntable 
rotations  with  respect  to  the  NUC  tower. 
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1 1. -2  Experimt'nt  Procedure 

The  sinRlc-axis  fluxgate  ^radiometer  was  mounted  on  a turntable  with 

the  center  of  rotation  for  both  the  angles  0 and  ip  located  at  the  center  of 

the  gradiometer.  The  three  axis  fluxgate  magnetometer  was  also  mounted  on 

the  turntable  near  the  center  of  rotation.  The  instruments  thus  measure 

the  magnetic  field  and  one  gradient  component  effectively  at  the  pivot  point 

of  the  turntalile.  In  addition,  one  output  of  the  fluxgate  gradiometer  provided 

a single  field  component  measurement,  and  could  be  used  as  a check  on  one  axis 

of  the  magnetometer.  The  entire  turntable  mechanism  could  then  be  attached 

to  any  one  of  several  posit ii>ns  on  a fiberglass  boom  which  is  located  on  the 

west  side  of  the  NUC  tower.  The  boom  may  be  positioned  at  different  elevations 

and  at  different  orientations  in  the  East-West  vertical  plane.  The  turntable 

mount  was  leveled  to  within  about  1 degree  during  measurements  by  comparison 

with  the  horizon.  This  procedure  was  checked  with  a transit.  Ttie  x^  and  x^ 

i j-oordinates  of  the  pivot  point  of  the  turntable  could  be  accurately  determined 

'f 

■ to  about  1 or  2 cm,  the  Xj  coordinate  could  be  estimated  to  about  10  cm, 

although  the  error  was  probably  much  less  than  this.  Coordinates  of  the  pivot 
point  for  the  nine  measurement  locations  (indicated  in  figure  1)  are  given 
In  table  1. 

1 The  turntable  could  be  remotely  oriented  (i.e.,  from  the  tower  deck)  at 

various  tilt  angles,  0,  by  turning  an  aluminum  pole  which  was  rigidly  attached 
' to  the  turntable  mechanism  and  supported  along  the  length  of  the  fiberglass 

boom  by  several  sleeves.  The  tilt  angle  could  be  precisely  positioned  at 
values  of  0,  i 22  1/2”,  and  ± 45”  by  aligning  a brass  pin  attached  to  the 
pole  near  the  turntable  with  one  of  five  similar  pins  located  on  the  sleeve 
nearest  the  turntable.  The  error  in  this  alignment  was  probably  no  more 
than  l”.  For  a fixed  tilt  angle,  full  36o”  rotations  in  the  angle  ip  could 
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be  accomplished  from  the  tower  deck  by  a pulley  and  cord  arrangement. 

The  turntable  could  be  locked  at  rotation  angles  of  every  15*^,  with  an 
uncertainty  of  less  than  1*^. 

Data  collection,  at  each  measurement  position  indicated  in  figure  1, 
proceeded  as  follows.  With  the  turntable  fixed  at  a tilt  angle  of  0°,  data 
from  the  gradiometer  and  each  axis  of  the  magnetometer  were  recorded  on  a data 
sheet  for  each  rotation  angle  of  30°.  This  data  was  collected  for  the  full 
360°  range,  including  both  the  0°  and  (redundant)  360°  positions,  and  then 
returned  to  the  0°  rotation  angle  and  the  data  recorded  for  this  postion 
again.  This  provided  data  at  12  distinct  rotation  angles,  with  the  0° 
rotation  angle  measured  3 different  times,  providing  a check  on  the  repeati- 
bllity  of  the  measurement  procedure.  With  the  turntable  returned  to  the  0° 
rotation  angle,  the  assembly  was  turned  to  the  + (or  - ) 22  1/2°  tilt  angle 
position  and  the  data  from  the  outputs  recorded  for  0°  rotation  only. 

The  turntable  was  then  turned  to  the  + (or-)  45°  tilt  position.  The  data 
from  each  instrument  was  recorded  at  every  30°  rotation  angle  for  the  0° 
to  360°  range.  In  the  same  manner  as  the  0°  tilt  measurement,  including  the 
repeated  measurements  at  0°  rotation  angle.  With  the  rotation  angle  fixed 
at  0°,  the  turntable  was  returned  to  the  + (or  - ) 22  1/2°,  and  then  to  0°  tilt 
position,  and  the  data  recorded  for  those  angles,  providing  additional  checks 
on  the  repeat Ibillty.  This  procedure  was  then  repeated  for  values  of  tilt 
angles  opposite  in  polarity.  This  entire  measurement  sequence  provided  data 
at  38  physically  distinct  angular  orientations  with  multiple  measurements 
at  each  tilt  angle  (for  the  rotation  angle  fixed  at  0°)  and  required  approxi- 
mately 1 hour  to  complete  for  each  boom  position  Indicated  In  figure  1 (after 
the  boom  was  In  position). 
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Point 

Label 

X2 (meters) 

(meters) 

A- 3 

-22.20 

-25.13 

A- 2 

-22.43 

-23.06 

A-1 

-22.66 

-20.92 

B-3 

-19.20 

-25.39 

B-2 

-19.43 

-22.92 

B-1 

-19.66 

-20.71 

C-3 

-16.20 

-25.11 

C-2 

-16.43 

-23.14 

C-1 

-16.66 

-20.79 

Table  1.  Coordinates  of  the  points  at  which  magnetic  field  and 
field  gradient  data  were  obtained.  (The  x^-coordlnate 
In  each  case  Is  0.) 
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III . Data  Analysis  and  Results  for  the  Magnetic  Field  and  Field  Gradients 
III-l  Instrument  Models 

The  output  of  a perfect  gradlometer,  as  a function  of  the  angles  0 
and  ipt  would  be  of  the  form  given  by  equation  (9).  For  a real  gradlometer, 
deviations  from  the  functional  form  of  (9)  can  be  expected.  For  small 
differences  in  the  areas  or  orientations  of  the  gradlometer  pick  up  loops, 
an  Imbalance  will  result  between  the  two  fluxgate  magnetometers  comprising 
the  gradlometer.  Additional  terms  in  the  output,  characteristic  of  a 
magnetometer,  will  then  be  present.  In  general,  the  output  of  the  single- 
axis gradlometer  is  a scalar  number,  and  may  be  written  as  a sum  of  scalar 
products.  Specifically: 


22 


6' 


•B  + x^'G'x^ 


(10) 


where  0^2  represents  any  dc  offset  in  the  output,  6* is  a vector  describing 
the  Imbalance  of  the  gradlometer,  and  the  last  term  is  the  output  of  a perfect 
gradlometer.  Higher  order  terms,  such  as  those  which  could  result  if  the 
magnetic  fields  varied  significantly  over  the  dimensions  of  a single  pick-up 
loop  in  the  gradlometer,  are  assumed  to  be  small. 

Equation  (10)  essentially  defines  a model  of  the  gradlometer,  which  is 
justified  only  if  the  data  can  be  adequately  described  in  terms  of  it.  It 
assumes  that  the  output  of  the  gradlometer  can  be  described  in  terms  of  the 
magnetic  field  B,  the  field  gradient  G and  four  parameters  characteristic 
of  the  instrument:  the  dc  offset  and  the  three  components  of  the  Imbalance 

vector  6. 

An  Instrument  model  may  also  be  developed  along  similar  lines  for  each 
axis  of  the  fluxgate  magnetometer.  Effects  which  were  included  in  the  analysis 
of  the  magnetometer  data  essentially  correspond  to  the  first  two  terms  appearing 
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on  the  right  hand  side  of  (10);  namely  the  output  of  one  magnetometer  component 
was  assumed  to  be  of  the  form; 

= C^'  + e(i)'  *B  , (11) 

where  i ” 1,  2,  or  3.  Three  instrument  parameters,  a dc  offset  and  two  com- 
ponents of  the  unit  vector  S(i)' , are  required  for  each  axis.  The  angular  depen- 
dence of  the  rotation  of  the  magnetometers  is  contained  in  the  term  S(i)'  *B, 
which  may  be  used  to  evaluate  the  components  of  the  magnetic  field  in  the 
fixed  coordinate  system.  For  rotations  described  by  an  operator  R,  the  magnetic 
field  transforms  as  RB.  Consequently,  for  the  specific  rotations  defined  by 
equations  (6)  through  (8),  the  term  £(i)'  *B,  may  be  written  as: 

t(i)'  *B  “ ej^(i)'((Bj  cos  6 - B^  sin  0)  cos  - B2  sin 

+ £2^^)'  t cos  0 - B^  sin  0)  sin  ip  + B2  cos 
+ e^d)'  [Bj^  sin  0 + B^  cos  0]  (12) 

For  an  axis  of  the  magnetometer  which  is  perfectly  aligned  in  the  primed 
coordinate  system  (i.e.,  is  precisely  aligned  with  respect  to  the  axial  gradi- 
ometer)  then  £(1) ' will  only  have  anx^'  component.  In  general,  the  alignment 
of  the  magnetometer  with  respect  to  the  gradlometer  was  quite  good,  and  the 
output  of  the  magnetometer  could  be  adequately  described  using  only  the  x^' 
component  for  S(l)' . However,  previous  work  had  shown  that  one  axis  was  off 
in  alignment  by  a few  degre<^s,  consequently  small  contributions  from  other 
components  were  Included  in  the  analysis  of  this  axis.  The  right  hand  side 
of  equation  (12),  with  the  £(i)'  components  replaced  by  the  5'  components  of 

the  imbalance  vector,  also  provides  the  functional  form  for  the  Imbalance 

term,  6**B,  appearing  in  equation  (10). 
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I1I-2  Determination  of  Field  Quantities 


In  order  to  determine  the  five  independent  elements  of  the  gradient 
matrix  from  the  data  collected  as  described  in  section  II-2,  a Fourier 
series  in  the  rotation  angle  was  used.  The  output  of  the  gradlometer 
was  expanded  as : 

N 

G'22  “ 2'  ^ (nij^)  + B^(6)  sin  (ni| ) ] . (13) 

n=l 

The  2N  + 1 coefficients  appearing  in  (13)  were  readily  calculated  using 
standard  Fourier  series  analysis.  Coefficients  up  to  and  Including  N*2 
are  expected  from  the  angular  dependence  predicted  by  equation  (10). 
Equating  the  coefficients  obtained  from  tie  data  with  those  arising  from 
the  various  terms  in  equation  (10)  allows  one  to  write  a set  of  equations 
for  the  elements  of  the  gradient  matrix  and  the  Fourier  coefficients  at 
the  three  tilt  angles  0 and  t 45°.  The  solution  to  these  equations  yields 
the  following  expressions  for  the  elements: 

Gii  - - + A2(-45°)]  - 2A2(0°)  . (14) 


G22  - I 1^2(45°)  + A2(-45°)1  . 


Gi2  4 62(0)  +-^(82(45°)  +82 (-45°)] 


Gi3  - A2(-45°)  - A2(45°) 


<=23  • ^ - B2(45»)I 
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The  result  of  this  analysis  yields  the  gradient  elements,  at  each  of 
the  positions  measured,  given  In  table  2.  Two  sets  of  gradient  elements 
are  given  for  the  position  labelled  A-1,  which  were  measured  on  different 
days  with  a complete  re-posltlonlng  of  the  boom  and  mount.  This  gives  a 
measure  of  the  repeatability  of  the  over-all  procedure  and  an  estimate  of 
the  uncertainty  In  the  results  presented.  A more  complete  analysis  of  the 
possible  errors  Is  contained  In  section  III-3. 

A similar  approach  was  used  for  the  calculation  of  the  magnetic  field 
components  from  the  magnetometer  data.  Again  the  output  of  a particular 
axis,  » was  expanded  as  a Fourier  series: 

1 ” 

®l'”2  ^ (ni(;)+b^(9)  sin  (nvp)  ] . (19) 

n»l 

Equating  the  coefficients  of  this  expansion  with  those  obtained  from 
equations  (11)  and  (12),  a set  of  equations  relating  the  magnetic  field 
components  to  the  Fourier  coefficients  is  obtained.  A set  of  equations 
analogous  to  (14)-(18)  results  for  each  axis  of  the  magnetometer.  Involving 
3q,  aj^,  and  bj^  at  the  various  tilt  angles,  6. 

The  results  of  this  analysis  yielded  values  for  all  components  of  the 
magnetic  field  from  each  axis  of  the  magnetometer  (with  one  exception,  which 
yields  only  two  components  of  the  field),  as  well  as  from  the  single  axis 
magnetometer  which  comprised  part  of  the  gradiometer.  In  addition,  certain 
axes  yield  more  than  one  value  for  a given  component.  All  totalled,  this 
analysis  gives  seven  distinct  measurements  of  the  component,  nine  measure- 
ments of  the  B2  component,  and  four  measurements  of  the  B^  component.  These 
independent  determinations  were  then  averaged  and  the  results  for  each  of 
the  nine  positions  where  data  was  collected  are  given  In  Table  3. 


19 


It  should  be  pointed  out  that  the  identification  of  the  Fourier  Series 
(13)  and  (19),  with  terms  appearing  in  the  instrument  model  equations  (10) 


j 


and  (11),  result  in  more  equations  than  unknowns.  Consequently,  the 
"solutions"  are  not  unique.  Additional  equations,  relating  the  Fourier 
coefficients,  field  quantities  and  instrument  parameters  are  also  obtained. 
The  equations  used  here  were  chosen  because  they  require  the  least  informa- 
tion concerning  the  instrument  parameters.  The  remaining  equations  not  used 
here  are  used  in  the  next  section  to  determine  the  instrument  parameters. 
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Point 

Label 


Si S2 S3 S2  S3 


A- 3 

-121.14 

-51.89 

-20.04 

335.26 

-47.24 

A- 2 

-111.16 

-58.61 

- 9.62 

320.93 

-94.22 

A-1 

-107,51 

-48.91 

-22.24 

286.48 

-132.87 

-100.96 

-46.38 

-18.53 

297.12 

-141.53 

B-3 

-226.02 

-84.13 

- 8,21 

527.96 

- 30.22 

B-2 

-224.81 

-80.96 

-15.63 

495.57 

-124.89 

B-1 

-176.94 

-82.93 

- 4.31 

432.80 

-210.85 

C-3 

-378.54 

-164.03 

-45.38 

875.56 

- 29.65 

C-2 

-333.02 

-141.13 

-22.92 

823.10 

-204.25 

C-1 

-301.13 

-143.24 

-34.31 

689.56 

-318.31 

Table  2: 

Magnetic  field  gradients  in  units  of  gamma/meter 

as 

measured 

at  nine  points  near 

the  NUC  tower.  Two 

measurements 

of  the 

gradients  were  made 

at  the  A-1  position 

on  March 

29 

(first  set) 

and  March  26  (second  set). 
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Point 

Label 

"l 

®2 

«3 

iBl 

A- 3 

22,278 

7,906 

41,150 

47,456 

A- 2 

22,223 

8,677 

40,327 

46,855 

A-1 

22,167 

7,239 

40,630 

46,847 

21,928 

7,195 

40,384 

46,513 

B-3 

21,374 

8,892 

41,675 

47,673 

B-2 

21,808 

8,221 

40,800 

46,987 

B-1 

21,237 

8,313 

40,504 

46,483 

C-3 

21,636 

11,032 

41,280 

47,894 

C-2 

21,877 

10,823 

40,588 

47,362 

C-1 

22,094 

9,450 

39,713 

46,417 

Table  3. 

Magnetic  field  components  and 

total 

magnetic 

field 

in  units 

of  gaimna  as 

measured  at  nine 

points 

near  the 

NUC  tower . Two 

measurements 

were  made  at  the 

A-1  position  on 

1 March 

29  (first 

set)  and  March  26  (second  set). 
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in- 3 Determination  of  Instrument  Parameters  and  Quality  of  Data  Fits 
The  additional  equations,  relating  the  Fourier  coefficients  from 
(13)  and  (19)  to  a combination  of  field  quantities  and  instrument  parameters, 
were  used  to  calculate  values  for  the  instrument  parameters.  For  the 
gradiometer  parameters,  the  following  relations  are  obtained: 


"22 


2-/2 


Ao<45°)  4 A„(-45“)  - /r  A„(0)  + ^ 


(20) 


^^T  B 


Aq(45°)  - Aq(-45°)  + G^3 


«3  = 


or 


Aq(0)  - 4G33  - «5G22  - C22 


(21) 


A^(0)Bj^  - 8^(0)  B2 

2 2 

«1  «2 


or 

A,  (45°)  (B,  - B.)  - B.(45°)B_  . 

1 ^13  1 2 


1s(B3  - 83)^  + 82^ 

or 

A^(-45°)  — (Bj^  + B3)  - 8^ (-45°)  B2 


4(B^  + 83)^  + B2^ 


(22) 
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A^(0)B2  + 


2 2 
Bi  +B^ 


or 


6! 


A (45°)B,  + 6.(45°)  — (B,  - B_) 
1 2 1 ^13 


4(Bj^  -B^)^  + 82^ 


A^(-45°)B2  ^i(-^5°)  — (B^  + B3) 



ij(B^  + B^)^  + 82^ 


(23) 


With  the  magnetic  field  quantities  determined,  the  instrument  para- 
meters could  then  be  calculated  from  any  Independent  set  of  4 of  these 
equations.  The  remaining  equations  should  then  be  satisfied,  within  the 
experimental  error  of  the  data  contained  on  the  right-hand  side.  An  optimum 

set  of  parameters  (ard  field  quantities)  could  also  be  determined  based  on  a 

2 

best  fit  analysis  such  as  minimizing  the  X . For  the  purposes  here,  however, 
a simple  average  of  the  Instrument  parameters  determined  from  all  of  the 
equations  was  performed.  Consequently,  the  dc  offset,  ^22*  was  uniquely 
calculated  from  (20);  the  two  values  of  6^  arising  from  (21)  were  averaged; 
and  the  three  values  for  each  of  6^  ^nd  coming  frcm  (22)  and  (23)  were 
averaged.  This  yields  the  parameters  given  in  Table  4 for  each  of  the 
positions  near  the  NUC  Tower  where  data  was  collected. 

As  is  clear  from  the  table,  the  parameters  which  model  the  output 
of  the  axial  gradlometer  varied  considerably  over  the  course  of  the 
experiment.  The  precise  reason  for  this  is  not  known,  but  is  probably 
associated  with  thermal  drifts  in  sensor  heads  or  electronics.  The  dates 
on  which  the  measurements  were  made  for  each  point  are  indicated  in  the 
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Point  C'22  (gamnia/meter)  6^  (meter  6^  (meter  6'^  (meter 

Label 


,-3 

12.78 

.00004 

-.00989 

-.00085 

l-2 

36.45 

.00021 

-.01023 

-.00086 

i-1 

00 

-.00005 

-.01027 

-.00090 

19.50 

.00005 

-.02182 

-.00098 

B-3 

9.35 

.00018 

-.02473 

-.00065 

B-2 

0.33 

.00024 

-.02451 

-.00060 

B-1 

10.46 

.00014 

-.02479 

- . 00080 

C-3 

3.47 

-.00006 

-.00809 

-.00064 

C-2 

0.73 

.00024 

-.01324 

-.00114 

C-1 

6.05 

.00008 

-.01498 

-.00115 

Table  4.  Gradlometer  instrument  parameters  as  determined  at  each  of  the  nine 
points  near  the  NUC  Tower  where  data  was  collected.  Dates  of  the 
measurements  were;  24  March  C-3;  25  March  C-2,  C-1;  26  March  A-1 
(second  set);  29  March  A-3,  A-2,  A-1  (first  set);  30  March  B-3,  B-2, 
B-1. 
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table  caption.  No  systematic  effect  in  the  dc  offset  is  apparent,  although 


the  sign  is  always  the  same.  However,  whi  le  the  imbalance  vector  was  relatively 
constant  on  a given  day,  it  changed  significantly  from  one  day  to  the  next. 

The  size  of  this  effect  may  be  readily  estimated  from  the  change  in  the 
dominant  component,  6^.  it  varied  from  about  0.01  to  0.02  (meter  ^), 
indicating  that  apparent  gradients  of  the  order  of  500  to  1000  gama/raeter 
could  be  observed  if  the  gradlometer  were  aligned  parallel  to  the  earth's 
field.  The  importance  of  Including  these  Imbalance  effects  in  understanding 
the  output  of  the  gradlometer  is  also  illustrated  by  a comparison  of  the 
data  obtained  or  two  different  days  at  the  same  measurement  point. 

Figure  7 shows  the  raw  data,  and  the  fits  based  on  the  gradlometer  model 

described  in  Section  III-l,  for  the  data  collected  on  26  and  29  March  at  ^ 

the  A-1  position.  The  output  of  the  gradlometer  is  significantly  different,  1 

several  hundred  gamma  /meter  variation  at  the  same  angular  orientation, 

although  the  gradient  matrix  elements  extracted  from  the  analysis  (given  ^ 

in  Table  2)  differ  by  no  more  than  9 gamma/meter . 

The  standard  deviations  in  the  fits  to  the  raw  data  using  the  functional 
form  of  equation  (10)  were  calculated  for  each  of  the  nine  measurement  points. 

With  the  field  quantities  from  Tables  2 and  3,  and  the  instrument  parameters 
given  in  Table  4,  the  expected  gradlometer  output  based  on  equation  (10) 
was  calculated  for  each  angular  orientation  at  which  data  was  collected. 

Letting  0^2  represent  the  actual  recorded  output  for  an  orien- 

tation labelled  by  a,  and  0^2 (exp)^  represent  the  actual  recorded  output 
for  char  orientation,  then  the  standard  deviation  is  given  by: 

I N f 

” I iPT  ^ ®22 

( o - 1 ^ 
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Kigure  6.  Plots  of  the  raw  data  (points)  and  fits  (solid  curves)  for  the 
gradiometer  output  at  the  A-1  position.  Units  are  gamma/ me  ter. 
For  each  tilt  angle,  the  larger  peak-to-peak  amplitude  curve  is 
for  the  26  March  measurements,  the  smaller  is  for  29  March. 


0 tilt  angle 


-AS  tilt  angle 


All  38  physicrlly  distinct  orientations  at  which  data  was  recorded  were 
Included  in  calculating  the  standard  deviation  for  each  measurement  point. 
These  standard  deviations  are  given  in  Table  6 and  provide  an  indication 
of  the  uncertainty  in  the  gradient  elements  determined  at  each  point. 

The  instrument  parameters  for  the  magnetometer  were  determined 
in  a similar  manner.  The  dc  offset,  for  each  axis  is  given  by: 

C[  - Aq(45°)  + Aq(-45°)  - Aq(0°)  . (2; 

2-/2~ 

The  dc  offsets  determined  this  way  are  given  in  Table  5.  The  g (i)'  were 
assumed  to  be  unit  vectors  in  the  x|  direction,  except  for  1=1  where 
three  components  were  calculated.  These  were  average  values  based  on  the 
data  from  all  measurement  points.  The  following  values  were  used  in  the 


calculation  of  the  standard  deviation  to  the  magnetic  field  fit: 


e(l)' 


i ^ 1 ; 


.99845  x| 


.044  ^2  + -034  , i = 1. 


This  corresponds  to  a 3.2  misalignment  of  the  i = 1 axis  of  the  magnetometer. 

The  standard  deviation  to  the  magnetic  field  fit  was  calculated  by 
summing  the  squares  of  the  vectorial  difference  between  the  actual  measured 


output  ard  the  expected  output  from  the  Instrument  model  (11): 


1 

K ^ 5' 

a - 1 


(exp)^  - b'  (cal)jj  I I 


This  gives  an  overall  statistical  error  estimate  for  the  total  magnetic 
field,  each  component  generally  being  more  accurate  than  this  total  standard 


deviation.  The  results  for  are  given  in  Table  6. 
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The  values  for  O,,  and  o arc  used  in  the  next  section  in  the 
development  of  the  analytic  field  model  ol  the  N'.  (.  tower.  They  provide 
a measure  of  the  relat've  uncertainty  associated  with  t lie  analysis  of 
the  data  in  terms  of  the  instrument  itKidels.  However,  the  standard  devia- 
tions calculated  here  only  in<  ludc'd  statistical  fluctuations  in  the  data, 
possible  systc'matic  errcirs  such  as  arise  friin  uncerta  int  ies  in  calibration 
coefficients  have  not  been  includ«*d.  effects  may  rc'siilt  in  an  over- 

all 10-207  uncertainty  in  the  magnet ir  field  and  gradient  values  given  in 
Tables  2 and  3. 

The  values  for  given  in  Table  b Indicate  that  the  instrument 
model  for  the  gradiometer  developed  in  Section  Ill-I  provides  a good  basis 
for  interpreting  the  raw  data  collected  during  the-  course  of  the  experi- 
ment. Typical  uncertainties  in  the  gradient  components  at  the  15-20  gamma/ 
meter  level  can  be  expected.  The  repeat  measurement  at  the  A-1  position 
suggests  that  it  may  be  even  somewhat  lower.  The  errors  in  a single  data 
point,  which  may  be  estimated  from  the  multiple  recordings  at  certain 
angles  as  discussed  in  Section  II-2,  are  typically  5 gamma/meter.  The 
standard  deviations  In  the  fits  to  the  data  are  only  modestly  larger  than 
this.  Barring  possible  systematic  errors,  these  results  Indicate  that 
the  gradient  components  given  in  Tab'e  2 are  probably  accurate  to  within 
about  one  standard  deviation  as  given  In  Table  6. 

The  standard  deviations  for  the  magnetic  field  fits,  o , given 

D 

in  Table  6 suggest  that  uncertainties  In  the  total  field  measurements 
may  be  as  large  as  700-1100  gamrvi.  The  repeat  measurement  at  the  A-1 
position  suggests  that  actual  uncertainties  for  the  magnetic  field  are 
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I 


somewhat  smaller.  The  total  field  difference  at  the  A-1  position  u'^ea- 
sured  on  two  different  days  was  only  about  300  gamma.  This  is  the  typical 
variation  on  the  total  field  observed  from  the  multiple  measurements  of 


a single  data  point  as  well.  These  results  indicate  that  the  magnetic 
field  data  given  in  Table  3 are  probably  at  least  as  accurate  as  the 
standard  deviations  given  in  Table  6,  again  barring  systematic  error. 

This  discussion  of  the  possible  errors  associated  with  the  deter- 
mination of  the  various  magnetic  field  quantities  also  demonstrates  two 
important  aspects  of  this  type  of  measurement.  The  uncertainties  asso- 
ciated with  measuring  the  magnetic  field  variation  from  point-to-point, 

as  indicated  by  the  standard  deviations  a_,  are  considerably  larger  than 

o 

the  actual  magnetic  gradients  underlying  the  variation  from  point-to- 
point.  This  is  the  reason  that  a gradiometer,  which  directly  measures 
spatial  variation  in  the  field,  was  used;  rather  than  trying  to  extract  the 
variation  from  an  array  of  magnetometer  data  alone.  Secondly,  the  use  of 
high  statistics  allowed  a reliable  determination  of  all  the  field  quan- 
tities of  Interest.  The  collection  of  data  at  38  physically  distinct 
angular  orientations,  together  with  the  known  rotation  properties  of  the 
gradient  matrix,  allowed  a precise  analysis  of  the  data  which  provided 
good  results  for  all  five  Independent  gradient  elements.  If  data  had  been 
collected  at  only  five  angular  orientations,  which  in  principle  should 
allow  the  determination  of  the  five  Independent  gradients,  the  Instrument 
parameters  would  have  been  completely  Indeterminate  and  several  hundred 
gamma/meter  errors  would  have  appeared  in  the  results.  The  design  of  the 
experiment  to  overdetermine  all  the  field  quantities  of  Interest,  resulted 
in  much  more  reliable  results  than  would  otherwise  have  been  the  case. 
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r* 


Point 

Label 

"1 

"3 

"2 

ni 

N3  • 

A- 3 

1672 

1204 

-959 

201 

A- 2 

1836 

2375 

-939 

104 

A-1 

1682 

1677 

1400 
2 382 

-1003 

-1208 

429 

-1127 

B-3 

1891 

1106 

-1279 

-70 

B-2 

240f 

1125 

-1125 

-304 

B-1 

1655 

1116 

-619 

429 

C-l 

1847 

924 

-817 

157 

C-2 

1514 

331 

-1418 

-228 

C-  1 

1443 

885 

-711 

206 

Table 

5.  Magnetometer 

dc  offsets  in  units 

of  gamma  as 

determined  at 

each  of  the 

nine  points  near  the 

NUC  tower  where  data  was 

collected . 

is  the  dc  offset  c 

if  the  single 

axis  field 

output  of  the  gradiometer. 
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Point 

Label 

Oj,  (gamma/meter) 

Og  (gamma) 

A- 3 

15 

774 

A-2 

25 

916 

A-1 

15 

707 

18 

1447 

B-3 

13 

854 

B-2 

21 

1178 

B-’ 

14 

660 

C-3 

12 

646 

C-2 

36 

1051 

C-1 

17 

759 

Table  6.  Standard  deviations  for  the  fits  to  the  raw  data  using  the 
Instrument  models  described  In  the  text. 
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IV.  AnaJ^^lc  _Model  of_  Nl]C  Tov^r 
I V- I Model  Descr i pt Ion 

An  analytic  model  of  the  magnetic  field  near  the  NUC  tower  was 
developed  in  order  to  provide  a simple  method  for  estimating  the  magnetic 
field  and  field  gradients  at  points  in  and  near  the  region  where  measure- 
ments were  made.  Three  candidate  models  were  considered.  One  included  a 
uniform  (earth's)  field  contribution  and  a contribution  from  a magnetic 
monopole.  The  second  included  an  additional  monopole,  and  tht  third  an 
additional  dipole.  The  third  model,  consisting  of  the  earth's  field,  a 
monopole,  and  a dipole,  provided  the  best  overall  fit  to  the  data,  and  this 
model  is  described  here.  The  following  definitions  are  used  in  this  section 


Bq  = earth's  magnetic  field,  independent  of  position.  (28) 


B„(r)  * M — , = magnetic  field  at  the  position  r,  of  a magnetic 

- I I J 

monopole  of  strength  M located  at  r^^.  (29) 


Bd(D 


3(r  - rj^)(r  - r^)  • D 


IE  - Ei)l 


l£  - Ed' 


magnetic  field  at  the  position  r,  of  a magnetic  dipole  of 


dipole  moment  D located  at  r^. 


(30) 


The  field  model 
matrix  elements,  G^j 


is  simply  the  sum  of  these  three  contributions, 
resulting  from  the  model  are: 


The 
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I 


'IJ 


TImI 


- 3 --- 


Ir-Ewi' 


3 I 

k-ol'  1 

IJ  " 1 ,2 

L J 

+ 


^r-rp>i  <r-r 


r-r. 


(31) 


The  four  parameters  of  the  monopole,  M and  r„,  and  the  six  parameters  of 

-M 

the  dipole,  D and  r^,  were  determined  by  a least  squares  fit  to  the  gradient 
elements  given  in  Table  2,  Each  point  was  weighted  according  to  the  standard 
deviations  given  in  Table  6.  Specifically, 


r-= 


G^.  (Equation  31 ) - G^^  (exp) 
Og  (exp) 


(32) 


was  summed  over  the  9 positions  and  5 independent  elements  and  a search 

2 

over  the  10  parameters  was  made  to  minimize  this  x • Tl>®  best  fit 

parameters  determined  in  this  manner  are  given  in  Table  7. 

With  the  monopole  and  dipole  contributions  to  the  field  model  determined, 

the  earth's  field  components  were  adjusted  to  obtain  a best  fit  to  the  mag- 

2 

netlc  field  data  given  in  Table  3.  Again  a X function  using  the  standard 
deviation  of  the  field  data  (o  from  Table  6)  was  minimized.  The  earth's 

D 

field  so  determined  Is: 


Bq  - [22,161  Xj  + 3,912  X2  + 41,491  x^J  gamma 


(33) 
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Monopole  Parameters 


1 7 2 

M * - ^ [2.14  X 10  I gamma-meter 

r„  - [5.13  X,  - 0.59  X - 26.09  x ] meters 


Dipole  Parameters 


1 Q 7 3 3 

D - 7—  [1.06  X 10  X,  + 3.13  X 10  X,  - 7.6  x 10  x,|  gamma-meter 

- 4tt  1 c 3 

tp  ” [-0.36  Xj  - 1.18  X2  - 23.90  x^l  meters 

Table  7.  Best  fit  parameters  for  the  field  model  using  a weighted 

least  squares  fit  to  the  measured  gradient  matrix  elements. 
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I V- 2 . Comparison  of  Model  with  Measured  Data 

The  m.ignetic  field  gradient  elements  calculated  at  each  of  the  nine 
measurement  positions  are  given  in  Table  8.  Below  each  entry  is  the 
difference  between  the  experimental  value  and  the  field  model  value. 

Table  9 gives  the  same  Information  for  the  magnetic  field  components. 

The  results  presented  in  Table  8 show  that  the  bulk  of  the  features 
exhibited  by  the  measured  data  are  reproduced  by  the  field  model  des- 
cribed here.  Both  the  monopole  and  dipole  term  play  an  important  role 
in  the  model.  For  example,  all  of  the  signs  of  the  off-diagonal  gradient 
elements  cannot  be  reproduced  by  a monopole  model  alone.  If  the  dipole 
moment  is  set  to  zero,  the  signs  of  all  the  Gj^2  el<?n'Pnts  are  positive, 
in  clear  contradiction  with  the  measured  data.  The  differences  between 
the  measured  data  and  the  values  calculated  from  the  model  are  not 
always  within  the  one  standard  deviation  uncertainty.  Gradient  differ- 
ences of  the  order  of  3o  are  apparent.  This  is  not  too  surprising, 

SJ 

however.  The  NUC  Tower  is  a very  complicated  structure;  simple  models 
of  the  field  could  only  be  expected  to  reproduce  most  of  the  gross  features. 
However,  most  of  the  larger  discrepancies  are  associated  with  the  larger 
elements,  so  that  as  a percentage  the  errors  in  the  model  predictions 
are  not  expected  to  be  very  large. 

The  results  of  the  model  for  the  magnetic  field  components  given 
in  Table  9 also  reproduces  all  of  the  bulk  features.  This  provides  a 
self-consistency  check  on  the  field  model  and  provides  confidence  that 
no  large  systematic  effects  have  been  overlooked. 
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The  overall  (uncorrelated)  X for  the  gradient  matrix  fit  is  85 

(45  data  points,  10  parameter  fit).  For  the  magnetic  field  fit,  the 
2 

X is  10  (15  data  points,  3 parameter  fit).  These  are  respectable  fits. 

A comparison  with  previous  data  taken  at  the  NUC  tower  and  also  with  other 
information  on  the  earth's  nvignetic  field  near  San  Diego  is  presented 
in  the  Appendix. 
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Point 

Label 


l X III. 

Label 

Si 

S2 

S3 

S2 

S3 

A- 3 

-161.15 

-27.03 

-12.34 

350.09 

-11. 35 

(40.01) 

(-24.86) 

(-7.70) 

(-14.83) 

(-35.89) 

A- 2 

-153,98 

-30.02 

-8.32 

326.71 

-62.14 

(42.82) 

(-28.59) 

(-1.30) 

(-5.78) 

(-32.08) 

A-1 

-143.22 

- 30 . 40 

-3.83 

287.07 

-102.66 

(35.71) 

(-18.51) 

(-18,41) 

(-0.59) 

(-30.21) 

B-3 

-241,00 

-60,38 

-26.88 

538.53 

-3.95 

(14.98) 

(-23.75) 

(18.67) 

(-10.57) 

(-26.27) 

B-2 

-228. 32 

-68,97 

-14.85 

492.37 

-113.70 

(3.51) 

(-11.99) 

(-0,78) 

(3.20) 

(-11.19) 

B-1 

-208.25 

-67.15 

-3.39 

412.84 

-181.92 

(31.31) 

(-15.78) 

(-1.15) 

(19.96) 

(-28.93) 

C-3 

-377.52 

-163.89 

-61,74 

881.71 

-23.41 

(-1.02) 

(-0.13) 

(16.36) 

(-6.15) 

(-6.24) 

C-2 

-357.15 

-177.00 

-31.73 

800.17 

-196.52 

(24.13) 

(35.87) 

(8.81) 

(22.93) 

(-7.73) 

C-1 

-318,41 

-165.40 

2,57 

629.77 

-331.47 

(17.28) 

(22.16) 

(-36.88) 

(59.79) 

(13.16) 

Table  8. 

Gradient 

matrix  elements 

In  gamroa/meter 

as  determined 

1 from  the 

best  fit  model  of  the  NUC  tower's  magnetic  field,  and  differences 
between  the  experimental  values  and  model  values. 
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Point 

Label 

|b| 

A- 3 



22,037 



7,740 

41,386 

47,522 

(2A1) 

(166) 

(-236) 

(-66) 

A-2 

22,022 

7,584 

41,018 

47,170 

(201) 

(1.093) 

(-691) 

(-315) 

A-1 

22,016 

7,334 

40,696 

46,847 

(131) 

(-95) 

(-66) 

(0) 

B-3 

21,918 

9,049 

41,422 

47,729 

(-3A4) 

(-157) 

(253) 

(-56) 

B-2 

21,881 

8,779 

40,733 

47,064 

(-73) 

(-558) 

(67) 

(-77) 

B-1 

21 ,877 

8,342 

40,247 

46,562 

(-640) 

(-29) 

(257) 

(-79) 

C-3 

21,601 

11 ,114 

41,277 

47,895 

(35) 

(-82) 

(3) 

(-1) 

C-2 

21,547 

10,696 

40,359 

46,984 

(330) 

(127) 

(229  ) 

(378) 

C-1 

21,555 

9,892 

39,527 

46,096 

(539) 

(-442) 

(18  6) 

(321) 

Table 

9.  Magnetic  field 

components  and  total 

field  In  units  of 

gamma  from 

the  best  fit  model  of  the  NUC  tower's  magnetic  field,  and  the 


difference  between  the  experimental  values  and  model  values. 
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Equivalent  Dipoles  of  Measured  Gradients 

Because  a steady  magnetic  field  In  free  space  Is  both  nondlvergent 
and  Irrotatlonal,  Its  gradients  at  a point  In  free  space  are  equivalent 
to  the  gradient  field  of  a magnetic  dipole  located  on  a sphere  of  unit 
radius  about  the  point.  As  we  describe  elsewhere^,  we  find  the  orientation 
and  location  of  an  equivalent  dipole  by  requiring  eigenvalues  and  principal 
axes  of  Its  gradient  field  to  be  respectively  equal  to  and  coincident  with 
eigenvalues  and  principal  values  of  gradients  at  the  field  point.  Because 
the  sense  of  an  eigenvector  defining  direction  along  a principal  axis  Is 
Indetermlnant , the  four  sets  of  eigenvectors  (e^^,  e^,  e^) , ®2’  » 

~^2’  ^®1*  ~®2’  define  the  same  principal  axes  and  so 

correspond  to  the  same  set  of  eigenvalues  X2*  X^).  As  a result,  each 

of  four  dipoles  located  on  a sphere  of  unit  radius  about  a field  point  give 


the  same  gradients  at  the  field  point. 


Figure  7 shows  orientation  and  location  of  a dipole  m^  that  produces  a 
gradient  field  whose  eigenvalues  are  Xj^,  X2»  and  X^»  and  whose  principal  axes 
coincide  with  the  set  of  eigenvectors  (a^,  @2*  ^3^*  dipole  axis  and 

position  vector  r^^  lie  In  the  plane  normal  to  the  eigenvector  position  ®2- 
Dashed  lines  In  Figure  7 delineate  orientation  and  location  of  dipoles  that 
produce  gradient  fields,  having  eigenvalues  X^^,  X2»  and  X^,  %rt>ose  principal 
axes  coincide  with  the  sets  of  eigenvectors  (8^^,  ~®2’  ~®3^ ’ ^”^1'  ®2*  ~®3^ ’ 


and  “®2*  ®3^  noted.  As  Is  evident,  the  four  equivalent  dipoles  lie 

in  the  plane  normal  to  the  eigenvector  ^2*  and  they  form  the  two  pairs 


and  (52,  £2;  -52. 


tan  X 


-r2) . The  relation 
0 < X ^ ■'f 


AO 


■1 


determines  Che  angle  x between  the  dipole  axis  and  position  vector  of 
an  equivalent  dipole,  and  the  relation 


U 3m 
o 

, 4 
4iTr 


1^2  + 


gives  magnitude  of  its  dipole  moment. 

Table  10  lists  the  eigenvalues*  and  polar  angles  defining  directions 
of  the  set  of  eigenvectors  (e^^,  e^,  e^)  for  each  measurement  position. 

Figure  8 shows  directions  of  the  set  of  eigenvectors  (e^^,  e^, 
frame  of  reference  {x^}  used  to  define  measured  gradient  elements.  The 
principal  axis  corresponding  to  eigenvalue  Is  approximately  aligned 
with  the  horizontal  x^  the  principal  axis  corresponding  to  eigen- 

value A^,  with  the  horizontal  x^^  axis;  and  the  principal  axis  corresponding 
to  eigenvalue  A^,  with  the  vertical  x^  axis.  The  plane  containing  axes 
of  equivalent  dipoles  is  nearly  vertical  with  x^  axis.  The  plane  con- 
taining axes  of  equivalent  dipoles  is  nearly  vertical  with  its  normal 
vector  e^  facing  north. 

Table  11  lists  dipole  strengths  and  polar  angles  defining  the  directions 
of  the  position  vector  r^^.and  of  the  dipole  axis  of  the  equivalent  dipole 
m^  located  on  a sphere  of  unit  radius  about  each  measurement  position. 

Figure  9 depicts  locations  and  orientation  of  the  equivalent  dipole  m^ 
in  the  vertical  plane  containing  measurement  positions.  Lengths  of 
dipole  vectors  are  drawn  in  proportion  to  their  strength. 


* Eigenvalues  are  indexed  so  that  A^  > A^  > A^. 
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Position 


Position  ffl/r  , amp/m 


VI.  Conclusion 

The  experiment  described  In  this  report  provided  a precise  measure- 
ment of  the  magnetic  field  and  associated  magnetic  field  gradients  at 
selected  points  In  a region  to  the  west  of  the  NUC  Oceanographic  Research 
Tower.  The  results  of  those  measurements  provide  a reliable  data  basis 
on  which  to  examine  Interference  and  noise  effects  for  the  proposed 
Internal  Wave  Magnetic  Sensing  Experiment.  An  analytic  model  of  the  NUC 
Tower's  magnetic  field  was  constructed  which  reproduces  the  Important 
features  of  the  data.  This  model  may  be  used  for  estimating  the  magnetic 
field  and  associated  gradients  In,  and  near,  the  region  In  which  data 
was  collected. 
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APPENDIX  A 

Comparison  with  Previous  Measurements 

A previous  survey  of  the  magnetic  field  near  the  NUC  tower.  In  a 
region  which  partially  overlaps  that  examined  by  this  report,  was  con- 
ducted In  January  1975  and  Is  described  In  another  report  (Physical 
Dynamics,  Inc.  report  PD-SM-76-112,  Appendix  A,  1976).  Data  were 
collected  In  that  experiment  using  a three  axis  fluxgate  magnetometer 
similar  to.  that  utilized  In  the  work  described  In  this  report.  No 
gradlometer  or  angular  rotation  data  were  collected  In  the  earlier  sur- 
vey. The  results  of  those  measurements,  together  with  the  field  values 
predicted  by  the  analytic  model  of  NUC  tower's  magnetic  field  described 
In  Section  IV  of  this  report,  are  given  In  Table  A-1. 

The  results  presented  In  the  table  Indicate  that  the  overall  features 
of  the  data  are  approximately  described  by  the  analytic  field  model  (e.g., 
the  trend  of  Increasing  In  magnitude  as  x^  becomes  less  negative) 
although  detailed  agreement  Is  lacking.  Some  of  the  disagreement  between 
the  and  B^  components  could  be  explained  by  a misalignment  of  the  magneto- 
meter with  respect  to  the  X2  axis  (axis  of  the  boom),  but  not  all  of  the 
discrepancy  can  be  described  by  this  effect  (e.g.,  the  total  field  does 
not  agree  at  each  point) . As  might  be  expected,  the  largest  discrepancies 
In  the  total  field  occur  at  the  two  points  which  are  the  farthest  from  the 
region  where  the  analytic  model  was  fit,  those  points  Indicated  as  near 
E-1  and  E-A. 

The  total  earth's  magnetic  field  extracted  from  the  data  taken  in 
this  experiment  (Eq.  33)  Is  about  47,200  gamma.  This  is  somewhat  smaller 
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than  the  value  of  about  49,000  ganma  given  by  Knecht  [David  J.  Knecht, 

"The  Geomagnetic  Field  (A  Revision  of  Chapter  11,  Handbook  of  Geophysics 
and  Space  Environments) Air  Force  Cambridge  Research  Laboratories 
AFCRL-72-0570,  1972,  pages  46-50]  for  the  San  Diego,  California  area 
(adjusted  to  1976) . The  vertical  component  given  In  Eq.  33  Is  very  close 
to  Knecht *8  value  of  41,500  gamma,  but  the  horizontal  component  (22,500 
gamma)  Is  smaller  by  several  thousand  gamma  (Knecht  gives  26,100  gamma, 
both  vertical  and  horizontal  components  adjusted  to  1976).  The  declination 
angle  for  the  field  of  Eq.  33  Is  only  about  10°  , compared  to  an  expected 
value  closer  to  14°.  However,  the  actual  alignment  of  the  coordinate 
system  defined  In  figure  5 with  respect  to  true  north  was  not  verified. 
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Coordinates  of 
Measurement  Point 


Magnetic  Field  - First  line  Is  experi- 
mental data;  Second 
line  Is  analytic  Model 
result 


X2  (meters)  x^ 

(meters) 

«1 

»2 

"3 

|b| 

- 9.9 

(~E-1) 

-19.6 

21,470 

17,258 

16,080 

14,986 

33,090 

33,023 

42,597 

40,161 

-12.9 

(~D-1) 

-19.6 

22,100 

20,399 

10,360 

12,067 

36,750 

37,001 

44,117 

43,941 

-15.9 

(~C-1) 

-19.6 

23,470 

21,430 

8,850 

9,907 

37,590 

38,887 

45,190 

45,492 

-21.9 

(-A-1) 

-19.6 

25,280 

21,987 

7,390 

7,393 

36,660 

40,419 

45,140 

46,602 

-21.5 

(~A-3) 

-24.2 

24,810 

22,003 

6,910 

7,973 

37,120 

41,184 

46,179 

47,369 

-11.5 

(~E-4) 

-28.1 

23,200 

20,686 

15,820 

16,340 

42,310 

45,015 

50,780 

47,889 

Table  A-1.  Magnetic  field  components  and  total  field 

In  units  of  gamma  as  measured  at  six  points 
near  the  NUC  In  January  1975,  and  the 
results  of  the  analytic  field  model.  The 
approximate  locations  In  terms  of  the  posi- 
tions Indicated  In  figure  1 are  shown  in 
parenthesis  (Position  E-4  Is  a few  meters 
above  position  E-3).  This  experimental  data 
is  courtesy  of  Richard  Zlemer  of  Physical 
Dynamics,  Inc.,  Santa  Monica,  CA. 
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METKIC  SYSTIM 


BASE  UNITS; 


Quantity 

Unit 

SI  Symbol 

Formula 

length 

metre 

m 

nMM 

kilogram 

ko 

time 

second 

s 

electric  current 

ampere 

A 

thermodynamic  temperature 

kelvin 

K 

amount  of  aubatance 

mole 

mol 

luminous  inteiuity 

candela 

cd 

SUPPLEMENTARY  UNITS; 

plane  angle 

radian 

rad 

solid  angle 

steradian 

sr 

DERIVED  UNITS; 

Acceleration 

metre  per  second  squared 

m/s 

activity  (of  a radioactive  source)  disintegration  per  second 

(disintegration  Vs 

angular  accelaralion 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

area 

square  metre 

m 

density 

kilogram  per  cubic  metre 

kg/m 

electric  capacitance 

farad 

K 

A-sA/ 

electrical  conductance 

siemens 

S 

AA/ 

electric  field  strength 

volt  per  metre 

V/m 

electric  inductance 

henry 

H 

V-s/A 

electric  potential  diffaienca 

volt 

V 

W/A 

electric  resistance 

ohm 

VIA 

electromotive  force 

volt 

V 

W/A 

energy 

joule 

1 

N-m 

entropy 

joule  per  kelvin 

1* 

force 

newton 

N 

kg-m/s 

frequency 

hertz 

Hz 

(cycle)/s 

illuminance 

lux 

lx 

Im/m 

luminance 

car.  dels  per  square  metre 

cd/m 

luminous  flux 

lumen 

Im 

cd-sr 

magnetic  field  strength 

ampere  per  metre 

A/m 

magnetic  flux 

weber 

Wb 

V-s 

magnetic  flux  density 

tesla 

T 

Wbhn 

magnetomotive  force 

ampere 

A 

power 

watt 

W 

V* 

pressure 

pascal 

Pa 

N/m 

quantity  of  electricity 

coulomb 

C 

A-s 

quantity  of  heat 

joule 

1 

N-m 

rwliant  intanaity 

watt  per  steradian 

W/sr 

specific  heat 

joule  par  kilogram-kelvin 

Vkg-IC 

stress 

pascal 

Pa 

Nhn 

thermal  conductivity 

wall  per  metre-kelvin 

W/m-K 

velocity 

metre  per  second 

... 

m/s 

viscosity,  dynamic 

pascal-second 

... 

Pa-s 

viscosity,  kinamatic 

square  metre  per  serxmd 

... 

m/s 

voltags 

vdl 

V 

W/A 

volume 

cubic  metre 

... 

m 

wavenumber 

reciprocal  metre 

... 

(wave)/m 

worfc 

joule 

1 

N-m 

SIPREraES: 

Multiplkatkm  Factors 

Prefix 

SI  Symbii 

1 000  000  000000  « in» 

ten 

T 

t 000  000000  - 10* 

mm 

C 

1 000  000  • 10* 

maga 

M 

1 000  • to* 

kilo 

k 

too  - 10< 

btsclo* 

h 

10  - 10’ 

dabs* 

da 

0.1  - to-' 

dad* 

d 

om  - io-> 

oatOi* 

c 

0.001  - to-* 

miUl 

m 

oiMnooi  • to-* 

■than 

«• 

0.000  000  001  - I0-* 

nano 

a 

OMt 

000000001  • 10*  •> 

pta» 

r 

OMtOOO 

OOOpOOOOl  - 10-’« 

ImMu 

f 

0000  000  000  OOO  000  oot  « to  ** 

aho 

a 

